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SUMMARY 

A linear finite strip plate element based on Mi ndl i n/Rei ssner plate theory 
is developed. The analysis is suitable for both thin and thick plates. In the 
formulation new transverse shear strains are introduced and assumed constant in 
each two-node linear strip. The element stiffness matrix is explicitly formu- 
lated for efficient computation and computer implementation. Numerical results 
showing the efficiency and predictive capability of the element for the analy- 
sis of plates are presented for different support and loading conditions and a 
wide range of thicknesses. No sign of shear locking phenomenon was observed 
with the newly developed element. 


1. INTRODUCTION 

The finite strip method, a combination of finite element and Fourier- 
series expansion, has important advantages for the analysis of a wide range of 
plate problems. The first finite strip bending element was developed following 
classical Kirchhoff thin plate theory (refs. 1 and 2). This formulation is 
very reliable for thin plate analysis but lack the ability to account for shear 
deformation which can be important when a plate becomes thick. In addition, 
displacement based finite element implementation of Kirchhoff plate theory 
requires ^-continuity (ref. 17, p. 172); this leads to high order or noncon- 
forming elements which is generally undesirable. 

Another theory, that has been subject of most recent research in finite 
plate elements, is the Mind! i n/Rei ssner theory (refs. 3 and 4). In the Mindlin 
theory, transverse shear deformation is included. Therefore Mindlin plate 
theory is applicable to modeling of classical thin plates as well as moderately 
thick, sandwich and composite plates (refs. 5 to 8) . The basic assumption is 
that a straight line originally normal to the middle surface of the plate 
remains straight but not necessary normal to the middle surface during deforma- 
tion. The vertical displacement is assumed not to vary through the plate 
thickness. Consequently, rotations are treated as independent variables. The 
finite element implementation requires only C 0 -continuity (ref. 17, p. 172). 
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Despite its mathematical elegance, overstiff numerical results are obtained 
when using linear and quadratic elements for the analysis of thin plates. This 
problem is very common and often refers to as a "shear locking" phenomenon. 

Several techniques have been proposed to remove the shear locking effect 
(refs. 9 and 10). However in the literature of finite strip, only two tech- 
niques are widely used. First, the use of "reduced integration" schemes for 
evaluation of the transverse shear components of the- element stiffness matrix 
is applied by Mawenya and Davies (ref. 6) for three-node quadratic Mindlin 
strip element. Benson and Hinton (ref. 11) used the same technique but 
extended the analysis to vibration and stability problems. The second tech- 
nique is the uses of polynomials of sufficiently high order to approximate the 
displacements. Onate and Suarez (ref. 12) used both techniques and compared 
the linear, quadratic, and cubic Mindlin strip elements using full, selective, 

and reduced integrations. Their conclusion is the two-node linear element with 

selective/reduced integration for the simplest and most efficient finite strip 
element. This agreed very well with the results of Hinton and Zienkiewicz 
(ref. 13). 

A new, simple two-node linear finite strip plate bending element for the 
analysis of very thin to thick structures is presented here. The essential 
ingredient of this new element lies in an assumed strain distribution technique 
introduced by McNeal (ref. 14) and Dvorkin (ref. 15) for developing finite ele- 
ment plate and shell elements. A new transverse shear strain is proposed and 
then constrained to equal the conventional transverse shear strain at selected 
poi nts . 

In section 2, the standard Mindlin finite strip plate bending element for- 
mulation is presented. The new formulation is given in section 3. Particular 

attention is focussed on points of evaluation for transverse shear strain. In 
addition to the theoretical development, numerical results for different sup- 
ports and loading conditions and a wide range of thickness are illustrated to 
assess the convergence and accuracy of the element are given in section 4. 


NOMENCLATURE 

A cross sectional of plate 

4* h 

strain displacement matrix for bending at node i for the a rn 
harmonic term 

strain displacement matrix for shear at node i for the 2.^ 
harmonic term 

Young 1 s Modulus 
k shear correction factor 

Ni shape function for node i 

q transverse loading per unit area 

n l l 

{R.} externally applied load vector at node i for the 2. tn harmonic 
1 term 



2 


t 

(uj) 

w 

»? 

3x 

I3y 

Yxz » Yyz 

e b 

e s 


natural coordinate variable 
thickness of plate 

nodal displacement vector at node i for the 2. th harmonic term 
transverse displacement 

nodal displacement at node i for 2. harmonic term 

section rotation in x-direction 

section rotation in y-direction 

transverse shear strains 

bending strain 

transverse shear strain 


for the fcth harmonic 
for the harmonic 

ah bending stress 

a s shear stress 

v Poisson's ratio 

ir total potential energy 


9 nodal rotation in x-direction at node i 

x i term 

% 

9 nodal rotation in y-direction at node i 

y i term 

k vector of curvature 


2. FORMULATION OF MINDLIN/REISSNER STRIP ELEMENT 

The finite strip method offers an alternative to the standard finite ele- 
ment method which uses polynomial functions in both x- and y-directions. 

The combination of polynomial and harmonic functions satisfy a priori the 
boundary conditions at the end of the strips and have the advantage of greatly 
reducing the number of equations to be solved for static analysis. For each 
type of supported condition at the strip ends, the necessary harmonic function 
is different. In this paper only simply supported conditions are considered. 

Consider a typical finite strip plate element, whose mid-surface lies in 
the x-y plane, with thickness t, longitudinal width b, strip width a, as 
shown in figure 1. The z-axis is normal to the plate mid-surface and the 
plate is assumed to be loaded by a pressure q. The sign conventions for posi- 
tive direction of the independent variables are illustrated in figure 1. 

For a two-node linear strip element, both the transverse displacement w 
and the section rotations |3 X and |3 y within an element are interpolated in 
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terms of nodal displacements w. and nodal rotations 0 , 6 as 
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r is the natural coordinate in x-direction 


and 



( 1 ) 


( 2 ) 


By using equation (1) simply supported boundary conditions are automati- 
cally satisfied at y = 0 and b. 
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3w 

ax 


ax ~ ay . ■ u 


Bending and shear strains are introduced separately as follow: 
Bending Strains (linear through the thickness) 


where 
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and 


is the strain-displacement matrix for bending at node i for the 


& harmonic term which can be written as 
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and {up is the nodal displacement vector as 


<“i> - *1 ®x, e 5 t J T 


Transverse Shear Strains (constant through the thickness) 
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3w 

ax + p x 


.ay ' P yJ SU1 1=1 
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where |^B.J is the strain-displacement matrix for shear at node i for the 

1.L S 

sr n harmonic term which can be written as 


ax S a 


s HirN.C^ 


- N 1 C <1 


The state of stress in the plate corresponds to a plane stress assumption, 
Considering an isotropic linear elastic material, the stress-strain relations 
are: 

[a.. 1 fl v 0 ] 


a b = a y " z 
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and 
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S 



_ E 

" 2(1 + v) e s 


where E is Young's modulus and v is Poisson's ratio. 

The total potential energy can be obtained as (ref. 17) 
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( 8 ) 


where 

t is the thickness of plate 
A is the cross-sectional area of plate 
q is transverse loading per unit area 

k is shear correction factor to account for nonuniformity of the transverse 
shear stresses through the plate thickness (usually set equal to 5/6) 

Substituting equation (1) and equations (4) to (7) into equation (8), the total 
potential energy can be rewritten as 


V = o 


k^D^ic dA + 




where 
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D - Etk 
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The loads are also resolved into a sine series in the y-direction similar 
to the displacement expression such that 
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< I 

q a s s 


(1 la) 


a=i 

where for a distributed load per unit area from y = c to y = d 

,d 


q a - b 


d(x,y)S^dy 


(lib) 


The strip equilibrium equations are obtained by imposing the stationary 
condition on ir, where w, |3 X , and |3y are independent variables. Using the 
orthogonality properties of the harmonic series (ref. 8) as 


EE KjK>-EE<> 

&=i i,j a=i i 


where 


KJ- KU« h M K]>JKI « 


for a single strip with node i and j 


< R t> - - ! 


q i[ N 


1 0 0 N 2 0 or dx 


’] T 


for node i under distributed load 
or 

{R-} = P sin 0p) 

for node i under a concentrated load P at y = c 


( 12 ) 
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The total resulting stiffness matrix is 


’k'l 

u 1 " 


R 1 

[k 2 ] 

u 2 


R 2 

k"l 

u n 


R n 


where 

[K ft ] is the global stiffness matrix for the 2.^ harmonic term 

{U a } is the global nodal displacement vector for the Q. th harmonic term 
and 

{R 2 '} is the global externally applied load vector for the harmonic term 

Equation (12) reflects the uncoupling of the stiffness matrix that leads 
to independent sets of equilibrium equations for each harmonic term. This fact 
can be taken advantage in simply supported case. That is, for other boundary 
conditions, the different harmonic terms may not always uncouple and a full 
stiffness matrix may have to be evaluated. 


3. ASSUMED STRAIN DISTRIBUTIONS 


The finite strip formulation presented in the previous section has some 
drawback. The element locks when the thickness is very thin using two-point 
Gauss quadrature integration for both bending and shear stiffness. This is due 
to the fact that the shear stiffness terms are overwhelming the bending stiff- 
ness terms and this leads to an overstiff element even though a very fine mesh 
is used. In the past this excessive stiffness may be attributed to "spurious 
shear" effects which were suppressed by reducing the order of the integration 
rule (ref. 9). To circumvent this locking phenomenon even though selective 
and reduced integration is a well-established approach, an assumed strain dis- 
tribution technique (refs. 14 and 15) is employed here. 

For each strip element, instead of using equation (5), new transverse 
shear strains are introduced as 


9. 1 / A C\ 

Y xz 2 vxz + Y xzJ 

2. 1 f B D \ 

Y yz 2 Vyz + Y yzJ 


( 17 ) 


where y‘ 


xz 


and 


Y , are the transverse shear strains evaluated at points 

A Z 
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and C with r = 0, y = b/4, and y = 3b/4 respectively (see fig. 2) using 
equation (5) as 


' xz 


xz 
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(18) 


with = sin(3Q/n74) , = si n(2.ir/4) , and J = a/2, and y!L . y V 7 are also 

Xj X/ y ^ y 

the transverse shear strains evaluated at points B and C with r at Gauss 
points and y = b/4, respectively (see fig. 2) using equation (5) as 


B r B 
Y yz C 2. 
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with cl = C° = cos(ftir/4). 


Hence, equation (17) can be rewritten with the new shear-strain displace- 
ment matrix as 
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Notice that the assumed shear strains are constant throughout the cross- 
sections and constrained to equal the shear strains of equation (5) at selected 
points. The choice of these points is of paramount importance in evaluating the 
predictive capability of the element even though the assumed strain distribu- 
tions in equation (17) are the integral part of the overall performance. With 
the shear strain displacement matrix of equation <20b) replacing equation (5b) 
the element stiffness matrix is formed explicitly as shown in figure 3 using 
exact order of integration. 


4. EVALUATION OF THE NEW FINITE STRIP ELEMENT 

The new finite strip element has been implemented into the finite element 

computer program FEAP (ref. 18) with relative ease. The subroutines written 
in FORTRAN 77 for the formulation of the element stiffness matrix consists of 

approximately 200 lines. All the results presented in this paper were obtained 

using double precision arithmetic on an IBM/PC-AT. 

Numerical problems have been tested in this section to show the perform- 
ance of the new element. The results are presented for different aspects: 
mesh size and harmonic term convergence characteristics, shear locking phenome- 
non as the thickness decreases, and shear force and bending moment prediction. 
Exact order of numerical integration is employed. Spurious zero energy mode 
effect (ref. 17) were not observed for the problems tested and the boundary 
conditions specified. Some of these results are compared with the analytical 
results obtained by other investigators (refs. 19 and 20). 


Shear Locking Investigation 

Using eight strip elements with four nonzero harmonic terms and v = 0.3, 
a simply supported square plate under two loading conditions: uniform and con- 
centrated loads, is investigated with a wide range of span/thickness ratio 
from 5 to 10^. The deflection at the center and quarter of the plate are plot- 
ted for both the full and reduced integrations of the troublesome transverse 
shear strain element, in figure 4 (uniform load case) and in figure 5 (concen- 
trated load case). Shear locking was not detected for the entire range of 
aspect ratios for the new element and reduced integration element. If the ele- 
ment was to exhibit shear locking as seen in the case of full integration of 
troublesome transverse shear strain, the displacements would be significantly 
less and less as the thickness decreased, due to overly stiff element behavior. 
Note that for the uniform load case the even harmonic terms are always zero as 
well as for the central concentrated load case. If however, the concentrated 
load was not located at the center, all the harmonic terms must be included. 


Convergence of Mesh Size and Harmonic Term 

To assess the convergence characteristics of the new element, a simply 
supported, uniformly loaded, square plate with t/b ratio equal to 0.01 and 
v = 0.3 is considered. For mesh convergence, the central deflections, central 
bending moments in both the x- and y-directions are tabulated in table I for 
four different meshes as well as the exact solution (ref. 19). Bending moments 
shown are linearly extrapolated from the Gauss points to the node. Also in the 
table, the number of degrees of freedom for each mesh is shown. This gives a 
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more realistic view of the computational cost. Graphical representation is 
shown in figure 6 by plotting percentage of error of the central deflections 
and central bending moment in the x-direction with mesh size. Good solutions 
are reached with only four strips and single harmonic term for both displace- 
ment and moment. The central deflection and bending moment for a four-node 
plate finite element (ref. 15) are also plotted in figure 6 for comparison pur- 
poses. This shows how rapid the rate of convergence is with no sign of shear 
locking. 

For the convergence of harmonic term, eight strip element mesh is used 
and the numerical results are tabulated in table II with four nonzero harmonic 
terms. Fairly good convergence is again obtained. Errors are less than 
0.4 percent when the third nonzero harmonic term is specified for both quanti- 
ties (see fig. 7). However the number of nonzero harmonic terms required for 
the analysis of more complicated problems may be increased. Note that this 
element is a low-order two-node strip. Numerical computation is minimum but 
convergence rate is relatively high. 


Predictions of Shear Force and Bending Moment 

For low-order element the standard finite strip formulation based on 
Mindl in/Rei ssner plate theory is well recognized to predict accurate displace- 
ments and fairly good bending moments when the selective and reduced integra- 
tion technique or high-order element is used (ref. 12). However, the shear 
force predictions are poor (see figs. 8 to 12) and hardly found in the finite 
strip literature even though they can be important in designing structures such 
as bridges and slabs. Therefore, the shear force predictive capability of this 
new finite strip element is studied and presented here. Due to the limited 
availability of analytical solutions, only four cases for shear forces and one 
case for bending moment are compared. 

The uniformly loaded square plates involving a variety of support condi- 
tions in the x-direction are investigated. In order to capture the steep gra- 
dient of the dependent variables near a plate edge, a rather fine mesh is 
modeled in the analysis with the new strip element. The results for varia- 
tions of shear forces and bending moments across the mid-plate in the various 
cases are plotted with analytical solutions by Kant and Hinton (ref. 20) in 
figures 8 to 12. These analytical solutions based on Mindlin plate theory 
assume the transverse displacement and sectional rotations similar to standard 
finite strip element. The system of differential equations is used and then 
numerically integrated using the so-called "segmentation method." The authors 
in that paper claimed that the analytical results compared favorably well with 
the finite strip method. 

The results of shear force and bending moment shown in figures 8 to 12 
along with the pertinent data and boundary conditions are in good agreement 
near the center of the plates. In the regions away from the center, the dif- 
ferences start to magnify and are approximately average at 15 percent near the 
edge of mid-plate. However, the curves for both solutions seem to follow the 
same pattern. Note that our element is only a simple, two-node linear element 
and its predictive capabilities (see figs. 8 to 12) are shown to exceed previ- 
ously developed elements. 
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5. CONCLUDING REMARKS 


A two-node linear finite strip element based on Mi ndl i n/Rei ssner plate 
theory is presented for static analysis of plates. A new shear strain distri- 
bution is assumed and connected to the standard shear strains at selected 
points. These points are chosen to remove shear locking phenomenon without the 
need for "reduced integration" technique. Due to the uncoupling nature of the 
finite strip method, the element stiffness matrix can be explicitly formulated 
for efficient computation and computer implementation. 

Several numerical studies were performed to assess the numerical perform- 
ance of the aforementioned finite strip element. Based on the results 
obtained, the following characteristics of the element can be stated: 

1. Simple, reliable, and efficient in computations 

2. Good convergence characteristics both in mesh size and harmonic term 

3. No shear locking effect for very thin situations 

4. Fairly accurate moment and shear force predictions 

5. Applicable to both thin and thick structures 
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TABLE I. - CONVERGENCE STUDY OF MESH SIZE FOR A SIMPLY 
SUPPORTED SQUARE PLATE UNDER UNIFORM LOAD 
WITH FOUR NONZERO HARMONIC TERMS 
[b/t = 100 and v = O.3.] 


No. of 
DOF 

No . of 
stri ps 

Central 
di spl acement , 
W( qb/D) 

Central bendi 

ng moments a 

Mx(qb) 

. My(qb) 

5 

2 

0.00348 

0.03691 

0.03865 

11 

4 

.00401 


.04699 

17 

6 

.00404 

Xii 

.04753 

23 

8 

.00405 

.04794 

.04766 

Exact 


0.00406 

0.0479 

0.0479 


a Linearly extrapolated from the Gauss points. 


TABLE II. - CONVERGENCE STUDY OF HARMONIC TERMS FOR A 
SIMPLY SUPPORTED SQUARE PLATE UNDER UNIFORM LOAD 
[b/t = 100 and v = O.3.] 


Harmoni c 
term 

Central 
di spl acement, 
W ( qb/D ) 

Central bending moments a 

Mx(qb) 

My(qb) 

a - 1 

-0.004101 

-0.04930 

-0.05162 

a = 3 

.000051 

.00158 

.00461 

a = 5 

-.000004 

-.00032 

-.00103 

a - 1 

.000001 

.00011 

.00038 

Sum 

-.004054 

-.04794 

-.04766 

Exact 

-0.00406 

-0.0479 

-0.0479 


a Li nearly extrapolated from the Gauss points. 
















FIGURE 1. - VIEW OF SINGLE STRIP ELEMENT AND NOTA- 
TIONS USED. 


y 



FIGURE 2. - SELECTED POINTS FOR ASSUMED SHEAR STRAINS IN 
EACH ELEMENT. ONE STRIP ELEMENT WITH x = a, AND 

y = b. 
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CENTRAL W. D/qb* 1 CENTRAL W, D/qb 


A 


NEW FINITE STRIP ELEMENT 
FULL-INTEGRATION FINITE STRIP ELEMENT 
REDUCED INTEGRATION FINITE STRIP ELEMENT 



LOG 10 b/t 



FIGURE 4. - SPAN/THICKNESS RATIO STUDY OF SIMPLY SUPPORTED SQUARE PLATE UNDER UNIFORM LOAD 
WITH FOUR NONZERO HARMONIC TERMS. 


O NEW FINITE STRIP ELEMENT 

•A FULL- INTEGRATION FINITE STRIP ELEMENT 

-O REDUCED INTEGRATION FINITE STRIP ELEMENT 




LOG 10 b/t 


FIGURE 5. - SPAN/THICKNESS RATIO STUDY OF SIMPLY SUPPORTED SQUARE PLATE UNDER CONCENTRATED 
LOAD WITH FOUR NONZERO HARMONIC TERMS. 
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FIGURE 6. - MESH CONVERGENCE STUDY FOR 
SIMPLY SUPPORT SQUARE PLATE UNDER 
UNIFORM LOAD ( V = 0.3) WITH ONE HAR- 
MONIC TERM. 


FIGURE 7. - CONVERGENCE STUDY WITH NUMBER 
OF NONZERO HARMONIC TERMS FOR SIMPLY 
SUPPORT SQUARE PLATE UNDER UNIFORM LOAD 
(V = 0.3). 


— NEW FINITE STRIP ELEMENT 

— ANALYTICAL (KANT AND HINTON) 

— ANALYTICAL (KIRCHHOFF) 



FIGURE 8. - VARIATION OF SHEAR FORCE 
Qxz ALONG Y/b = 0.5 FOR SQUARE 
PLATE UNDER UNIFORM LOAD WITH FOUR 
NONZERO HARMONIC TERMS, 


FIGURE 9. - VARIATION OF SHEAR FORCE 
Qxz ALONG Y/b = 0.5 FOR SQUARE 
PLATE UNDER UNIFORM LOAD WITH FOUR 
NONZERO HARMONIC TERMS. 
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Qxz/qb 


NEW FINITE STRIP ELEMENT 



FIGURE 10. - VARIATION OF SHEAR FORCE 
Qxz ALONG y/b = 0.5 FOR SQUARE 
PLATE UNDER UNIFORM LOAD WITH FOUR 
NONZERO HARMONIC TERMS. 


FIGURE 11. - VARIATION OF SHEAR FORCE 
Qxz ALONG y/b = 0.5 FOR SQUARE 
PLATE UNDER UNIFORM LOAD WITH FOUR 
NONZERO HARMONIC TERMS. 



FIGURE 12. - VARIATION OF MOMENT Mx ALONG 
y/b = 0.5 FOR SIMPLY SUPPORTED SQUARE PLATE 
UNDER UNIFORM LOAD WITH FOUR NONZERO HAR- 
MONIC TERMS (V = 0.3). 
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